Are neutron stars with crystalline colour superconducting cores interesting for LIGO? 
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We estimate the maximal deformation that can be sustained by a rotating neutron star with a 
crystalline colour superconducting quark core. Our results suggest that current gravitational-wave 
data from LIGO have already reached the level where a detection would have been possible over a 
wide range of the poorly constrained QCD parameters. This leads to the non-trivial conclusion that 
compact objects do not contain maximally strained colour crystalline cores drawn from this range 
of parameter space. We discuss the uncertainties associated with our simple model and how it can 
be improved in the future. 



Introduction. — A spinning neutron star may be an 
interesting gravitational-wave source provided that it is 
deformed in some way. Such asymmetries, colloquially 
referred to as "mountains" , may arise in a number of dif- 
ferent ways. The elastic crust of the star may sustain 
shear stresses [H, 0] , the size of which are limited by the 
breaking strain of the crustal lattice. A magnetic star 
in which the magnetic field is misaligned visa-vi the ro- 
tation axis may also be deformed in an interesting way 
[3, Estimates of the possible level of neutron star 
asymmetry are of great current interest given the improv- 
ing upper limits obtained by LIGO [3| . The best current 
limit set by the detectors corresponds to a deformation of 
e < 7.1 X 10"'^ in PSR 2124-3358, a neutron star spinning 
at 202.8 Hz. This should be compared to the maximal 
mountain size predicted from theory. We have recently 
concluded that the largest crust mountain one should ex- 
pect on an isolated neutron star is e « 2.4 x 10~^ 0]. This 
estimate is obtained assuming that the breaking strain of 
the crust is (ibr = 10~^, which may be optimistic. The 
comparison shows that the sensitivity of the detectors 
have reached the level where the observations are begin- 
ning to confront the theory. Of course, we are still some 
way away from testing more conservative models. Recall, 
for example, that the crust breaking strain is usually as- 
sumed to lie in the range (Jbr — 10~^ — 10~^. If the true 
value is at the lower end of this range, then we would 
still be far away from a detection. It is also important 
to keep in mind that the theoretical estimates are based 
on a maximisation argument. It is not at all clear that 
there are physical avenues that lead to a star being de- 
formed at this level. As a final caveat, we note that there 
also exist spin-down upper limits on ellipticities, where 
conservation of energy is used to produce a moment of 
inertia and distance dependent bound on e. However, 
these are less robust than the direct upper limits, and 
will not be our focus here. 

Despite the various caveats, the existing observational 
data are very interesting. As shown by Owen Q, maxi- 
mally strained solid quark stars could radiate at a level 
where a detection would have been possible. It may also 
be the case that strong emission could occur from an 



elastic phase in the neutron star core. Core deformations 
are, in fact, likely to be more significant than crustal 
ones, since the high density region provides a larger con- 
tribution to the quadrupole moment. The problem has 
been the lack of "quantitative" models. This may have 
changed recently, with the suggestion of crystalline colour 
superconducting phases in the deep core (see 0, S, i 
and references therein). Our aim with this letter is to 
adapt our recently developed framework for estimating 
the maximal crustal mountains in a neutron star [5] to 
consider elastic cores. This leads to a simple model for a 
neutron star with a crystalline core of deconfined quarks 
and a normal hadron fluid envelope. We compare the 
attainable mountains within this model to the current 
limits set by LIGO and ask to what extent observations 
are already confronting QCD modelling. 

The elastic core model. — As discussed in 0, @| it 
may be energetically favourable for colour superconduct- 
ing quark matter to form a crystalline structure. The 
shear modulus can be estimated as ^] 
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It depends on the gap parameter A and the quark chem- 
ical potential /ic, for which the authors of Q suggest 
estimated ranges of 

350MeV < fic< SOOMeV , 5MeV < A < 25MeV . (2) 

There is naturally a significant uncertainty associated 
with these estimates. In order to keep the analysis sim- 
ple, we will consider the shear modulus to be constant in 
the core and dependent only on the combination (A/Zc)^, 
a parameter on which one may then be able to place 
constraints from observations. In reality, the various pa- 
rameters are expected to be (weakly) density dependent, 
but it does not seem relevant to try to account for this 
in a first study. 

We consider a simple model of a neutron star with an 
elastic core of deconfined quarks and a fluid exterior. The 
equation of state for the exterior is taken to be an n = 1 
polytrope, while the core is described by an incompress- 
ible fluid. This is a reasonable first approximation for a 
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compact star with a deconfined quark core. The partic- 
ular advantage of this simple core-mantle model is that 
the Newtonian hydrostatic equilibrium equations 



(3) 



where p is the pressure, $ the gravitational potential and 
p the density, can be solved analytically. We match the 
solutions for the two regions at a given transition den- 
sity pc and impose the constraints that the total mass 
of the star is M = IAMq and the radius is i? 10 
km. These constraints are obviously not necessary, but 
they restrict the available parameter space. This way 
one can construct a sequence of stellar models by vary- 
ing the transition density, which we take to be in the 
range 3 < PcJ Pn < 8 (where p^ is the nuclear saturation 
density). Two examples of the resulting density distribu- 
tions are shown in Figure [1] together with the radius of 
the solid core as a function of the transition density. 

In order to calculate the maximum ellipticity we extend 
the formalism developed in 0]. In essence, we perturb 
the background configuration assuming that the pertur- 
bations have a 1^22 f) angular dependence. Deviations 
from sphericity build up strain in the core. It breaks 
when a — a^r at some point. Here, a is the modu- 
lus of the strain tensor as defined in As we have 
already mentioned, the breaking strain abr is highly un- 
certain. To make progress we will assume that it lies in 
the usual range considered for neutron stars, i.e. we take 
10-5 < CTbr < 10-2. 

To maximise the ellipticity we thus assume the maxi- 
mum strain in the core and solve the hydrostatic equilib- 
rium equations, which now include elastic terms due to 
the deformations: 



where we have defined the stress tensor 
t„ 



P ( Va6 + VbCa - g'JahV"^ 
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with p the shear modulus, which we consider constant, 
and the components of the displacement vector for 
the crystalline phase. As boundary conditions we require 
regularity at the centre of the star and continuity of the 
components of the traction both at the surface and at the 
interface between the solid and the fluid, where we also 
impose that the deformed shape of the star maximise the 
strain in the core, as described in [3]. 

Having solved for the equilibrium shape of the per- 
turbed star we can calculate the ellipticity 
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which, for an observer on the z axis, is tied to the 
gravitational-wave strain amplitude by 

^ = —J z — ' y^) 



where u is the neutron star spin frequency, Izz it's prin- 
cipal moment of inertia and r the distance from Earth 
Pi . In order to facilitate a direction comparison with the 
ellipticities discussed in [S], we fix I^z = lO^^g cm^ rather 
than calculating it for each of our stellar models. This is 
just a convention, since it is clear from ^ and ([7]) that 
Izz actually does not affect the final expression for the 
gravitational-wave amplitude. 

Results and discussion. — The ellipticities that we cal- 
culate depend on the breaking strain ahr and the com- 
bination (Apc)"^. The results also depend on the chosen 
value of the transition density between the core and the 
crust. This is obvious, as one would expect a more signifi- 
cant elastic core to be able to sustain larger deformations. 

Our main results are illustrated in the left panel of Fig- 
ure [21 where we show the obtained ellipticities for both 
the maximum breaking strain usually considered for a 
neutron star, a^r — 10~^, and a more conservative value, 
(Tbr = IQ-^. For each of these values we consider the 
region between the maximum and the minimum value of 
the shear modulus p for the range of parameters sug- 
gested in 0], cf. 0. We compare these ellipticities to 
the best current upper limit set by LIGO using the S3/S4 
science runs, which is e = 7.1 x 10"'^ for PSR J2124-3358 
0. It is quite clear that, if the maximum breaking strain 
applies and the star is maximally deformed, then LIGO 
would have made a detection. Such a detection would 
have implied, cf. Figure [21 that the transition to a quark 
core would have to take place below p « 7.5pn- The fact 
that no detection was made rules out such a scenario. Of 
course, this non-detection does not say whether it is the 
QCD model that is at fault (i.e. there is no such high 
density core) or whether the QCD core exists but simply 
isn't significantly strained. 

These results are obviously interesting, especially since 
the sensitivity of gravitational wave searches are set to 
improve. Looking ahead, the soon to be completed S5 
science run will provide approximately one year's worth 
of data at initial design sensitivity. The S3/S4 data 
spanned approximately one month with a noise floor at 
near twice design sensitivity, so an upper bound from the 
S5 run will give a maximum ellipticity a factor of 2 • \/T2 
tighter than for S3/S4 (the upper bound scales linearly 
with the noise floor and as the inverse square root of the 
observation duration). For pulsar J2 124-3358 this would 
correspond to a maximum ellipticity of approximately 
1 X IQ-^. This could allow detections even if the break- 
ing strain were at the lower end of the range we have 
considered. For instance, if a^r = IQ-^ then detection 
would be possible if the crust-core transition occurs at a 
density no greater than p « 5pn, assuming parameters 
at the upper end of the range given in ([2|) . Looking fur- 
ther ahead, one year of data from Advanced LIGO would 
yield a further improvement of a factor of ten in sensi- 
tivity, corresponding to e = 1 x IQ-^ for J2124-3358. A 
detection would then be possible providing the crust-core 
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FIG. 1: The left panel shows the density profiles for two different values of the ratio pc/pn- One can see quite 
clearly that (for a fixed mass and radius of IAMq and 10 km, respectively) the core becomes significantly smaller 
as the transition density rises. This fact is also illustrated in the right panel, which shows the radius of the core 
as a function of the transition density. 



transition density lies below p « 3pn regardless of where 
the shear modulus and breaking strain parameters lie in 
their respective likely ranges. A non-detection (i.e. up- 
per bound) would demonstrate that the star in question 
does not contain such a maximally strained core with the 
appropriate transition density. 

Our analysis suggests neutron stars containing colour 
crystalhne cores are promising candidates for gravita- 
tional wave detection, and that in the event of no de- 
tection being made the corresponding upper hmits place 
non-trivial constrains on the deep interiors of compact 
objects, showing that either the strain in the core is small 
or that the high density equation of state does not sup- 
port such solid phases. Of course, the model we have 
considered is very simphstic, and needs to be improved 
in a number of ways if we want the results to be truly 
quantitative. Most important would be to improve our 
understanding of the crystalline colour superconductor. 
We have assumed that it behaves like ordinary elastic 
matter with a given shear modulus. This is a natural first 
assumption, but it may not be the complete picture. For 
example, we have not accounted for the "superfluid" na- 
ture of the core. Yet this could be an important omission. 
Perhaps a comparison between the crystalline colour su- 
perconductor and supersolid Helium 

0, El [II would 

help improve our understanding? This is an interesting 
possibility since supersolid Helium is amenable to labo- 
ratory experiments. 

Up to this point we have assumed that the nature of 
the phase transition at the interface between the quark 
core and the hadron fluid is such that there is no dis- 
continuity in the density profile. However, there could 
be a (potentially sizeable) discontinuity at the interface, 
see for example Fig. 4 in [l^]. To gain some insight into 
how this might affect our results, we have carried out a 
comparison of the ellipticities of three stars that can be 
easily modelled within our framework. Firstly we consid- 



ered a star consisting of the usual solid core surrounded 
by fluid, calculating I^x — lyy (and hence the ellipticity) 
by integrating over the whole star in the conventional 
manner. Then we took the same star but integrated over 
the core only. Finally we took a bare core with no fluid 
and with the same mass and radius as the core in the 
first two models. The results of these calculations (as a 
function of core density) are shown in the right panel of 
Figure [21 

Interestingly, there is little difference between the 
quadrupole of the bare core and the quadrupole inte- 
grated over only the core of the core plus fluid star. How- 
ever, the addition of the fluid envelope serves to decrease 
the ellipticity, the decrease being more severe for stars 
with higher core-fluid transition densities. This decrease 
can be understood in a simple intuitive way. Unlike the 
solid, the fluid cannot support shear strains, so its pertur- 
bation from spherical symmetry is sourced entirely by the 
perturbed gravitational potential. The fluid's outer sur- 
face will therefore be close to spherical. The fluid there- 
fore serves to 'fill in' the depressions of the solid core, 
resulting in a partial cancellation between the Ixx — lyy 
contributed by the core and fiuid. Clearly, this cancella- 
tion will be more substantial the closer the density of the 
core and fluid, the case of a smooth density transition 
and a bare core being extreme cases. 

This leads us to conjecture that in the case of a sub- 
stantial density step at the core-envelope interface, the 
cancellation would be far from complete, leading to an el- 
lipticity closer to the bare core values shown in the right 
panel of figure [H For a small core, this could lead to an 
order of magnitude increase in the ellipticities calculated 
assuming a smooth density transition, increasing further 
the importance of gravitational-wave observations in con- 
fronting QCD. A rigorous treatment of such a star should 
employ a realistic equation of state (such as that em- 
ployed in 13*1) which would naturally incorporate any 
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Po/P. Pc/P„ 

FIG. 2: In the left panel, the ellipticities obtained from our analysis are shown as a function of the core transition 
density Pc/pn, where pn is the nuclear saturation density. We consider two breaking strains, a maximum value 
(Tbr = 10~^ and a minimum of (Tbr ~ 10~^. For each of these values we examine the region between the maximum 
and the minimum shear modulus p, obtained for the range of parameters in cf. (|2]). The shaded region 
indicates the region permitted by LIGO observations and the horizontal line is the current LIGO upper limit 
of e = 7.1 X 10"'^ for PSR J2124-3358 If the breaking strain is close to the maximum the observations are 
already confronting the theory, showing that J2124-3358 at least does not contain such a maximally strained core. 
If the breaking strain is close to the minimum, however, no such conclusion can be drawn. In the right panel we 
compare the ellipticity obtained by considering the full star, only the core of the full star and just the naked core, 
with no fluid around it. We take the breaking strain to be (Jbr = 10~^. As one can see the results for the two 
cores do not differ significantly, while the ellipticity of the core plus fluid star is smaller, more so as the core size 
decreases at higher transition densities. 



density step. 

The model we have considered is obviously simpUs- 
tic in many other ways. In particular, one would want 
to relax the assumption that A/^c is constant. In or- 
der to make real improvements, one would (again) want 
to use a realistic equation of state. This would require 
the calculation to be carried out within General Relativ- 
ity. Then the determination of the background solution 
is straightforward, but the calculation of the mountain 
size would require implementing the General Relativistic 
theory of elasticity, see for example flij. To date, there 
have been no such calculations. Work in this direction 
should clearly be encouraged. 

Finally, we need to improve our understanding of the 
breaking strain. The range of values that we have used, 
10"^ < (Tbr < 10~^, is relevant for a crust consisting 
of normal matter. However, the physics of the core is 
very different from that of the crust. There is no reason 
to believe that the estimates on the breaking strain for 
the crust should be applicable to the core. The response 
of the crystalline quark matter to large stresses is also 
uncertain. Normal matter will be predominantly brittle 
and break into pieces when the temperature is sufficiently 
far below the melting temperature and will respond by 
plastic flow (up to some limit) otherwise. How an elastic 
quark core will respond is completely unknown. Yet, for 
our purposes it may not matter which scenario is realised 
as long as the timescale for plastic flow at a given strain 
is longer than the observation time. 
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